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6 PSEUDO-CHAOTIC COMMUNICATION METHOD 

7 EXPLOITING SYMBOLIC DYNAMICS 

8 This invention was made with government assistance provided by the Army 

9 Research Office (DAAG55-98-1-0269). The government has certain rights in this 
10 invention. 
11 

1 2 FIELD OF THE INVENTION 

13 The field of the invention is data communication. The invention is 

14 particularly applicable to ultra- wide bandwidth impulse-radio communication systems. 
15 

1 6 BACKGROUND OF THE INVENTION 

17 The continually increasing reliance on wireless forms of communication 

18 creates reliability and privacy problems. Data should be reliably transmitted from a 

19 transmitter to a receiver. In particular, the communication should be resistant to noise, 

20 interference, and possibly to interception by unintended parties. 

21 In the last few years there has been a rapidly growing interest in ultra- wide 

22 bandwidth (UWB) impulse radio (IR) communication systems. These systems make use 

23 of ultra-short duration pulses that yield ultra- wide bandwidth signals characterized by 

24 extremely low power spectral densities. UWB-IR systems are particularly promising for 

25 short-range wireless communications as they combine reduced complexity with low 

26 power consumption, low probability of detection (LPD), immunity to multipath fading, 

27 and multi-user capabilities. Current UWB-IR communication systems employ pseudo- 
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1 random noise (PN) coding for channelization purposes and pulse-position modulation 

2 (PPM) for encoding the binary information. 

3 Others have proposed aperiodic sequences of pulses in the context of 

4 chaos-based communication system. Additional work has relied upon the self- 

5 synchronizing properties of two chaotic systems. In such a system, data is modulated into 

6 pulse trains using variable time delays and is decodable by a coherent receiver having a 

7 chaotic generator matched to the generator used in the transmitter. Such system is known 

8 in the art as a Chaotic Pulse Position Modulation (CPPM) scheme. 

9 Such chaotic dynamical systems have been proposed to address the 

10 problem of communication privacy. Chaotic signals exhibit a broad continuous spectrum 

1 1 and have been studied in connection with spread- spectrum applications. The irregular 

12 nature of a chaotic signal makes it difficult to intercept and decode. In many instances a 

1 3 chaotic signal will be indistinguishable from noise and interference to receivers not 

14 having knowledge of the chaotic signal used for transmission. In the context of UWB 

1 5 systems the use of nonperiodic (chaotic) codes enhances the spread-spectrum 

16 characteristics of the system by removing the spectral features of the signal transmitted. 

17 This results in a lower probability of interception/detection (LPI/LPD) and possibly less 

18 interference towards other users. This makes the chaos-based communication systems 

19 attractive. 

20 There remains a need for improved chaotic coding/modulation methods to 

21 produce such attractive communication systems. It is an object of the invention to meet 

22 that need. 
23 
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1 SUMMARY OF THE INVENTION 

2 The invention is a pseudo-chaotic coding/modulation method. The coding 

3 method of the invention exploits symbolic dynamics of a chaotic map at the transmitter to 

4 encode data. This produces an encoding system that synthesizes the chaotic map based 

5 upon the data to be transmitted. 

6 In a preferred embodiment, pseudo-chaotic iterates are generated from a 



7 digital implementation of a Bernoulli shift map. The output of the shift map is translated 

8 by a mapping, preferably implemented by a digital signal processor, to allow transitions 

9 between states in a transmitted signal to differ, and the translated map is used to drive a 

10 modulator (for example PPM, FSK, PSK, QAM, etc.). In the specific case of pulse- 

1 1 position modulation (PPM) the translated map is used to modulate pulse train positions 

12 within a periodic synchronization frame. The preferred embodiment uses a shift register 



13 to implement an approximation of the Bernoulli shift map acting as a form of 

14 convolutional code with a number of states equal to the symbolic states defined on the 

15 chaotic map. A receiver may use fewer states and still decode the data signal, allowing 

1 6 receiver scalability. 

17 A preferred transmitter accepts digital data for coding. The digital data is 

1 8 allocated to symbolic states according to a chaotic map. The pseudo-chaotically coded 

19 data is converted to analog form and modulated into synchronization frames in a 

20 transmitted signal. 

21 Another preferred embodiment accepts digital data for coding. The digital 

22 data is encoded by applying a chaotic map having N states, where N=2 M , with M being 

23 the number of bits in the shift register implementing the Bernoulli shift map. The 

24 symbolic states are defined on the chaotic map according to a Markov partition and the 

25 sequence of the states constitutes the encoder output. The pseudo-chaotically coded data 

26 is converted to analog form and modulated to produce a modulated signal for 

27 transmission. 
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1 The preferred embodiments thus include features to produce useful coding, 

2 decoding, and modulation methods. Artisans will understand the important features may 

3 be applied to communication systems in different ways while still realizing advantages of 

4 the invention. 
5 

6 BRIEF DESCRIPTION OF THE DRAWINGS 

7 Other features, objects and advantages of the invention will be apparent by 

8 reference to the detailed description and the drawings, of which: 

9 FIG. 1 is a block diagram of a preferred pseudo-chaotic coding/modulation 

10 method of the invention; 

1 1 FIG. 2 illustrates a Bernoulli shift map with the definition of the symbolic 

12 dynamics in accordance with a preferred embodiment of the invention; 

13 FIG. 3(a) illustrates implementation of the Bernoulli shift map through the 

14 M-bit shift register and D/A converter of FIG. 1 ; 

15 FIG. 3(b) illustrates the effect on state quantization by the shift register of 

16 FIG. 3(a) with M=4; 

17 FIG. 4 illustrates an effect of a Gray/binary translation from the Bernoulli 

1 8 shift map to a tent map; 

19 FIG. 5 illustrates a preferred periodic synchronization frame for the 

20 modulation conducted by the transmitter of FIG. 1 when pulse-position modulation is 

21 used; 

22 FIG. 6 shows the tent map having the same definition of the symbolic 

23 dynamics as in FIG. 2 and illustrating the symbolic states definition according to a 

24 Markov partition with N=8 states; 

25 FIG. 7 illustrates the transition state diagram associated to the dynamics of 

26 the tent map with the state definition illustrated in FIG. 6; and 

27 FIG. 8 illustrates the trellis diagram associated with the transition state 

28 diagram of FIG. 7. 

4 





1 



2 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
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Referring now to FIG. 1, a preferred embodiment system constructed 



4 according to the invention includes a transmitter/encoder 10 that communicates with at 

5 least one receiver/decoder 12 over a channel 14, typically (but not necessarily) a wireless 

6 channel. The encoder 10 accepts a data stream and transmits data into the channel 14 

7 using a pseudo-chaotic encoding. In the preferred embodiment, coding is conducted 

8 primarily by a shift register 16 that implements an approximation of a Bernoulli shift and 

9 a digital signal processor 18. The digital signal processor translates dynamics of the 

10 Bernoulli shift register, while permitting simple transformations, e.g. a Gray/binary 

1 1 conversion, or the realization of more complex pseudo-chaotic maps, which might be 

12 used for spectral shaping purposes and/or to enhance maximum likelihood detection by 

13 the receiver/decoder 12. Output from the digital signal processor is converted by a D/A 

14 converter and used to modulate data in a modulator 22. In a preferred embodiment, 

15 pulse-position modulation is used in a periodic synchronization frame. In general, the 

16 symbolic dynamics approach enables realization of a Viterbi detector 24 in the receiver 

17 with a number N R of states lower than the levels in the transmitter. 
18 

19 

20 I, Encoding 

21 A. Shift Map 

22 In the preferred embodiment of FIG. 1, the shift register 16 implements a 

23 shift of an incoming bit stream. In considering the preferred embodiment and the 

24 implemented shift, symbolic dynamics theory will be helpful in aiding understanding of 

25 the invention. 

26 We first consider the shift map and resulting symbolic dynamics. £ 2 

27 indicates the symbol space of binary "0"s and 'T's, that is 2 2 = {(s 0 SiS 2 . . . ):s; = 0 or s x = 

28 1 }. Then, the shift map a: Z 2 — » £ 2 is defined in Equation (1) as: 
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that is the shift map simply "forgets" the first entry in a sequence, and shifts all other 

3 entries one position to the left. Incidentally, the shift map is perhaps the simplest example 

4 of chaotic dynamics and it possesses all the peculiar features of chaotic systems. 

5 An equivalent way of studying the shift map is to represent the state x of 

6 the corresponding discrete (one-dimensional) dynamical system as a binary expansion, as 

7 in Equation (2): 

GO 

x = OAM3- = Z 2 ~ Jb j < 2 ) 

9 where each of the bits bj is either a "0" or a "1", and* e [0,1]. Then, the effect of the 

10 shift map applied to the binary sequence {bj}"} =l is described by the so-called Bernoulli 

1 1 shift map of Equation (3): 



x k+l = 2x k mod 1 (3) 



14 

15 A graph representing the Bernoulli shift map is shown in FIG. 2. Referring to 

16 Equation (3), successive iterates of x are obtained by moving the separating point one 

17 position to the right (multiplication by 2) and setting to zero the first integer digit 

18 (modulo operation). Hence, digits which are initially far to the right of the separating 

19 point, and thus have only a very slight influence on the value of x eventually become the 

20 first fractional digit. In sum, a small change of the initial conditions eventually makes a 

21 large change in x k , confirming the sensitivity to initial conditions of the chaotic Bernoulli 

22 shift. 
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2 B. Symbolic Dynamics 

3 Symbolic dynamics may be defined as a "coarse-grained" description of the evolution of 

4 a dynamical system. The idea is to partition the state space and to associate a symbol to 

5 each partition. Then, a trajectory of the dynamical system can be analyzed as a symbolic 

6 sequence. In the case of the Bernoulli shift map, shown in FIG. 2, the state space is 

7 represented by the invariant interval I = [0,1]. As mentioned, the first step for 

8 characterizing the symbolic dynamics of a given dynamical system consists of 

9 introducing a proper partition of the state space. To ensure that the symbolic dynamics 

1 0 give rise to a topological Markov chain, a so-called Markov partition has to be selected. 

11 In the case of the Bernoulli shift map (with N=2), a Markov partition may be selected by 

12 simply splitting the interval I = [0, 1] with respect to the critical point c = 0.5 and, 

13 correspondingly, we define the two subintervals I 0 = [0, 0.5) and i! = [0.5, 1), as 

14 illustrated in FIG. 2. 

15 In order to obtain a symbolic description of the dynamics of the chaotic 

16 map under consideration, FIG. 2 associates the binary symbol "0" to the subinterval I 0 

17 and the symbol "1" to the subinterval Ii. Then, the evolution of the state of the Bernoulli 

18 map can be described in terms of a symbolic sequence S = {010010 ... }. 
19 

20 C. Shift Register and DSP 

21 The basic idea behind the symbolic dynamics encoder 26, which includes 

22 the shift register 16, DSP 18, and D/A converter 20 in FIG. 1, is that the Bernoulli shift 

23 process may be implemented by means of an infinite length shift register R where at each 

24 step the most significant bit (MSB) is discarded. The shift operation corresponds to a 

25 multiplication by a factor 2, while discarding the MSB at each step is equivalent to a 

26 modulo 1 operation. 

27 A practical shift register, of course, has finite length and on the other hand 

28 it is impossible to specify the initial conditions with infinite precision. Encoding binary 
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1 data exploiting the Bernoulli shift map requires consideration of a finite-length shift 

2 register R that is fed with a binary data stream c(k) to be transmitted. We assign the most 

3 recent bit of data the least significant bit (LSB) position in the shift register 16. At each 

4 step (or clock impulse) the new bit of data is copied into the LSB position of the shift 

5 register while the (old) MSB is discarded. The situation is depicted in FIG. 3(a). 

6 We assume that the binary data stream c(k) feeding the shift register is an 

7 i.i.d. (independently identically distributed) sequence. A randomization of a data stream 

8 may be obtained, if desirable, by use of a data compressor and/or a data scrambler 28, as 

9 shown in the preferred FIG. 1 embodiment, to pre-code the data stream prior to entry into 
10 the shift register 16. 

11 

12 1 . State Quantization 

13 In practice, due to the finite length of the shift register R, the dynamics of 

14 the Bernoulli shift can only be approximated, as the admissible states assume only 

1 5 discrete values. Namely, by considering a M-bit shift register, the generic state x can be 

16 expressed as: 

M 

x = 0.b l b 2 ...b M =Y J 2 ~ Jb j ( 4 ) 

18 to be compared with Equation (2), where bj and b M represent the MSB and the LSB, 

1 9 respectively. Obviously, the approximation can be made arbitrarily precise by increasing 

20 the length M of the shift register. The effect of the state quantization on the Bernoulli 

21 shift map is shown in FIG. 3(b), for the case M = 4. 
22 

23 D. Control of Chaos in the Transmitter 

24 The preferred embodiment transmitter implements a form of predictive 

25 control with respect to the symbolic dynamics of the Bernoulli shift, with the state 
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1 definition given by Equation (4). The symbolic dynamics of the Bernoulli shift with the 

2 Markov partition are determined solely by the successive values of the MSB in the shift 

3 register R. Referring to Equation (4), the first bit on the right of the separating point 

4 determines whether the iterate falls within I 0 or l\ (the remaining bits deciding only the 

5 relative position within Ij with j = 0, 1). The value of the MSB at step k coincides with 

6 the value of the LSB M steps before, at step (k - M). In turn, the LSB contains the 

7 current bit of data c(k). In this sense, the scheme is predictive in its nature. 

8 The present invention solves the control problem a priori, by direct 

9 "synthesis" of the pseudo-chaotic signal, starting from the binary data to be transmitted. 

10 From this perspective, the injection of a new bit of data in the LSB of the shift register 16 

1 1 may be interpreted as a perturbation of the state of the dynamical system, in order to 

12 make it follow the desired symbolic sequence. The strength of the perturbation is of the 

13 order of 2 _M , thus can be made arbitrarily small by increasing the number M of bits of the 

14 shift register. 
15 

16 E. Map Shaping/Translation 

17 The digital signal processor 18, in FIG. 1, may be realized by a Gray/binary 

1 8 converter. The purpose of the Gray /binary conversion is to translate the dynamics of the 

19 Bernoulli shift map, given by Equation (3), into those of the tent map, described by: 
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x k+l =1-2 x k -0.5 (5) 



21 

22 

23 as illustrated in FIG. 4. This is done in order to achieve a greater robustness of the system 

24 in presence of noise by avoiding (zero-order) discontinuities in the map. The tent map has 

25 very similar dynamics to the Bernoulli shift map; in particular, both admit the same 

26 constant invariant density. Also, their symbolic dynamics with respect to the Markov 

27 partition can be related to each other. 
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2 In the preferred embodiment, however, the digital signal processor 1 8 is 

3 capable of generating more complex chaotic maps. This may be useful for spectral 

4 shaping purposes and/or for enhancing the Maximum-Likelihood detection. 
5 

6 

7 F. Modulation to Produce a Transmitted Signal. 

8 The signal from the digital signal processor is transformed into an analog 

9 signal by the D/A converter 20, which is used to drive the modulator 22 . In the case of 

10 PPM a constant offset is then added to this analog signal to form the modulation signal 

1 1 m(t) and used to drive a Pulse Position Modulator. Please note that the following two 

12 sections assume the modulator 22 conducts PPM. 
13 

14 1. Synchronization Frame 

15 In the preferred embodiment, each pulse is allocated, according to the 

16 pseudo-chaotic modulation signal, within a periodic synchronization frame, as shown 

17 schematically in FIG. 5. This assumes the existence of a periodic reference (with period 

1 8 T f ) such that only one pulse train (or in the simplest case a single pulse) for each user 

19 (User A or User B in FIG. 5) is transmitted within each frame time T f , coinciding with 

20 the symbol period. In FIG. 5, each frame time includes a "guard" time interval, t g , 

21 (proportional to the constant offset) for avoiding overlapping between adjacent pulse 

22 trains. Each pulse can occur at any of TV = 2 M (discrete) time instants, where M is the 

23 number of bits of the shift register 16 and of the D/A converter 20. In FIG. 5, the timeslot 

24 corresponding to each level of the pulse-position modulation has been denoted by t s . Note 

25 that the value of t s is limited from below by the system time-base resolution. 

26 Consider the case of data encoding in accordance with the symbolic 

27 dynamics of the tent map. Consider the exemplary partition with respect to the critical 

28 point c = 0.5, and indicate with t c the corresponding time delay from the beginning of the 
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1 frame. Then, by indicating with t k the relative time (again referred to the beginning of the 

2 frame) at which the k-th pulse train begins, if t k < t c a "0" is being transmitted, while a 

3 "1" is being transmitted if t k > t c . The situation is shown in FIG. 5. The use of a 

4 synchronization frame enhances the robustness of coding in the presence of noise and 

5 spurious pulses by preventing error propagation phenomena. 

6 2. Multilevel PPM 

7 Modulation conducted in accordance with the invention may be considered 

8 a N-PPM multilevel modulation with input bit "coded" through the pseudo-chaotic map 

9 with rate 1/M, where M = log 2 (N), realizing a sort of random coded modulation. For 

10 every user bit, a number of M channel bits are transmitted and detected. N R indicates the 

1 1 number of levels at the receiver. According to the invention, this may differ from the 

12 number N of levels at the transmitter. In particular, the relation N R <N holds, allowing a 

13 certain degree of freedom in the design of the receiver 12, and admitting classes of 

14 receivers having different complexity levels which are still capable of decoding a signal 

1 5 modulated in accordance with the invention. 
16 

17 

18 II. Decoding 

19 The receiver 12 decodes the signal from the channel 14. The receiver 12 

20 includes a demodulator 30, a decoder 32 matched to the chaotic map, an output mapper 

21 35 and a decompression/descrambling block 34. 

22 A. Pulse Position Demodulation 

23 In FIG. 1, the ideal detection of the incoming pulses (affected by noise) is 

24 achieved by a pulse position demodulator (PPD) 30 consisting of a pulse correlator 

25 matched to the pulse shape and a decision circuit. In the simple case of a single pulse per 

26 frame, the PPD 30 may be realized by an Integrate and Dump (I&D) filter which 

27 estimates the position of the pulse within each frame time. We assume that a normalized 

28 signal d(t) ee [0, 1] is available as output of the PPD. 
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2 B. Threshold Discriminator 

3 Demodulated data must be decoded by a decoder 32. In the simplest case 

4 the demodulator output d(t) can be decoded by a simple threshold decoder. The decoder 

5 32 is followed by an output mapper 35 and a decompressor/descrambler 34. A decision 

6 threshold 36 in the decoder 32 should be set according to the partition corresponding to 

7 N R = 2, for example to the value c = 0.5, coinciding with the critical point of the tent map 

8 as shown in FIG. 6. 

9 1 . Coding for Noise Resistance 

10 With the threshold detection, it is reasonable to expect that most of the error 

1 1 events will originate from pulses corresponding to values of the pseudo-chaotic iterates 

12 close to the partition point of the tent map (c = 0.5), separating the symbol "0" from the 

13 symbol "1". To reduce this particular error event probability a noise "gap" may be created 

14 in the chaotic map around its partition point c. This can be obtained by a proper pre- 

15 encoding of the data, establishing forbidden sequences of bits. One possibility is to add a 

16 code that avoids series of consecutive zeros. This constraint is usually known as run- 

17 length limit constraint and denoted by RLL (0, k), where k indicates the maximum 

1 8 number of allowed consecutive zeros. 
19 

20 C. Viterbi Detection/Maximum-Likelihood Estimation 

2 1 The Viterbi detector 24 performs a Maximum-Likelihood (ML) estimation 

22 of the transmitted sequence. Other techniques for deriving a sub-optimal estimator for a 

23 chaotic process in additive white Gaussian noise (AWGN) may also be used. 
24 

25 1 . The Tent Map as a Markov Chain 

26 Consider again the Markov partition of the invariant interval I = [0, 1] of 

27 the tent map in N "states". The situation is illustrated in FIG. 6, for the case N = 8. From 

28 FIG. 6 , transitions between different states are governed by the dynamics of the map. In 
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1 particular, only certain transitions (in this case two) are allowed from each state. For 

2 example, referring to FIG. 6 , it is clear that the interval corresponding to the state 1 can 

3 only map to itself or to the state 2. More precisely, the transition taking place depends on 

4 which subinterval associated with each state, the generic iterate x k belongs to. These two 

5 possible transitions have been labeled with 0 and 1, respectively. 

6 Now, the Markov chain associated with the tent map (with the partition in 

7 FIG. 6) can be represented by means of an equivalent transition diagram, as illustrated in 

8 FIG. 7. The transition branches have been labeled according to the finer partition in 

9 subintervals, as discussed above and illustrated in FIG. 6. FIG. 7 also shows (within a 

10 box) the value of the output corresponding to each state, according to the encoding of the 

1 1 data. The trellis corresponding to the transition diagram of FIG. 7 is shown in FIG. 8 
12 

13 2 . Detector Scalability 

14 An interesting feature of the invention is the possibility of realizing the 

15 Viterbi detector 24 with a number N R of states lower than the states used by a transmitter 

16 14 to encode data. This is obtained simply by matching the Viterbi detector 24 at the 

17 receiver 12 to the map with a lower number of states. Given the transmitter number of 

18 states, N, there is a broad range of possibilities for decoding the signal produced by a 

19 transmitter implementing the invention. Receivers having Viterbi detectors with 

20 complexities N R = 2, 4, 8, . . . N, are possible. This scalability property enables 



21 receivers having different complexities and performance to decode the same transmitted 

22 signal. 
23 

24 

25 3. Output Mapper/Decoding Other than Gray Conversion 

26 When a generic chaotic map is implemented by the DSP 1 8 in the 

27 transmitter 10, the receiver 12 should include a further function after the detector 32 in 

28 order to exploit the scalability feature. This additional function reconstructs the 
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1 transmitted message given the estimated sequence of states provided by the detector 32 

2 itself. This function is an output mapper, that will be, in general, a finite state machine 

3 (FSM). For particular cases and for appropriate choices of the input/output labels of the 

4 encoding map like in the Bernoulli shift and in the tent map, this FSM is not necessary. 
5 

6 III. Multi-user Access 

7 Artisans will appreciate that the invention is amenable of multi-user access. 

8 For example, in the case of PPM multiple-access of the channel may be realized by 

9 assigning different pulse trains for each user and correspondingly different matched 

10 filters at the receiver side. The situation is illustrated schematically in FIG. 5, showing 

1 1 pulse trains associated to two different users, User A and User B. Other multiplexing 

12 techniques, like for example time-division multiplexing (TDM), may be also utilized for 

13 multiple-access purposes. 
14 

15 While various embodiments of the present invention have been shown and 

16 described, it should be understood that other modifications, substitutions and alternatives 

17 are apparent to one of ordinary skill in the art. Such modifications, substitutions and 

18 alternatives can be made without departing from the spirit and scope of the invention, 

19 which should be determined from the appended claims. 

20 Various features of the invention are set forth in the appended claims. 
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